ON THE MORDELL WEIL GROUP OF THE ELLIPTIC CURVE 

y 2 = x 3 + ri 



YASUTSUGU FUJITA AND TADAHISA NARA 

Abstract. We study an infinite family of Mordell curves (i.e. the elliptic curves 
in the form y 2 = x 3 + n, n 6 Z) over Q with three explicit integral points. We 
show that the points are independent in certain cases. We describe how to compute 
bounds of the canonical heights of the points. Using the result we show that any 
pair in the three points can always be a part of a basis of the free part of the 
Mordell- Weil group. 



o 



1. Introduction 



H 

r-| . Let E be an elliptic curve over a number field K. It is known that the set of rational 

points E(K) is a finitely generated abelian group by the Mordell- Weil theorem. If the 
absolute value of the discriminant of E is not large, we can practically use Cremona's 
program 'mwrank'. However there is no known algorithm which determines the struc- 
ture of E(K) even if K — Q. The difficulties come from the free part of the group. 
We are interested in the families of elliptic curves of which we can at least partially 
determine the structure of the Mordell-Weil group, that is, the families which have 
explicit points which can be in a system of generators of the Mordell-Weil group. 
In the paper [3] Duquesne considered an infinite family of elliptic curves in the form 
y 2 = x 3 — nx. He showed that the curves in the family have two explicit integral 
points which can always be in a system of generators. Recently, the first author and 
Terai ([!]) generalized Duquesne's theorem on generators and showed that the same is 
true for infinitely many binary forms n = n(k, I) in 7L\k, /]. In this paper we consider 
an infinite family of elliptic curves in the form of y 2 = x 3 + n with three explicit 
$_j ■ integral points. 

Let a, b be integers and 

(1.1) E af) : y 2 = x 3 + a 6 + 16b 6 
the elliptic curve over Q. We put 

(1.2) P 1 = (-a 2 , 46 3 ), P 2 = (2ab, a 3 + 4b 3 ), P 3 = (-2ab, a 3 - 4b 3 ). 

Then it is easy to see that they are in E a ^(Q). In this paper we prove the following 
theorem. 

Theorem 1.3. Assume that a,b are relatively prime integers with a,b>3 such that 
a 6 + 16& 6 is square-free, ab is odd and b is divisible by 3 but not by 9. Then the rank 
of the Mordell-Weil group E at t,(Q) is at least 3 and any pair of two points {Pi,Pj} 
{i = 1, 2, 3, i 7^ j) can always be in a system of generators of E at b( 
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Remark 1.4. If n is square-free and not equal to 1, the elliptic curve y 2 = x 3 + n has 
no rational torsion points by P, Theorem 5.3]. Therefore Pi, P2, P3 are non-torsion 
in the situation of Theorem 11.31 

We prove Theorem 11.31 along similar lines to Duquesne's ([3]). First we com- 
pute bounds of the canonical heights of Pj's (i = 1,2,3). For that purpose we use 
the decomposition of the canonical height into the sum of local heights. The non- 
archimedean part of local heights is computed by Silverman's algorithm. To compute 
the archimedean part of local heights, we use Tate's series. Then the archimedean part 
is in the form log |x(Pj)| + | Yl^Lo^ - " 1 I^C^Pi)!) where z(P) is a certain function. 
But to do this, it is necessary that the ^-coordinates of the points on the curve are 
away from zero. To deal with the difficulty, we shift the elliptic curve in the direction 
of x-axis. We set the shifting width by functions of a, b such as 2a 2 + 4b 2 , 3a 2 + 4b 2 . 
Then the computation of the bounds comes down to a maximum- minimum problem 
of elementary functions. Further we find in our case z(P) above is bounded indepen- 
dently of a, b and P. Thanks to this, we obtain an upper bound and a lower bound 
whose difference is a constant. 

On the other hand we find a uniform lower bound of the canonical height. It is a 
bound independent of the points on the elliptic curve. This is computed by Cohen's 
algorithm. 

Next using those bounds we estimate the lattice indices of {p, P,} {i, j — 1, 2, 3, i ^ 
j) (for the definition of the lattice index see Section The key theorem is Siksek's 
theorem which comes from the theory of quadratic forms. 

The goal of the proof of Theorem 11.31 is to show that the lattice indices equal 1. 
By the estimete using Siksek's theorem we can show that the lattice indices are less 
than 5. On the other hand we prove that the lattice indices are not divisible by 2, 3 
by an argument of the descent. This completes the proof. 

The organization of this paper is as follows. In Section [2] we review basic notations 
of elliptic curves. We also review the canonical height and the local height function. 
In Section [3] we compute bounds of the canonical heights of Pi, P2, P3. In Section H] 
we compute an uniform lower bound of the canonical height. In Section Owe estimate 
the lattice indices by applying Siksek's theorem to the results of Sections |3] and HI In 
Section [6] we prove that the lattice indices are not divisible by 2, 3 by an argument 
of the descent. Then we complete the proof of Theorem 11.31 Further we prove that 
the family of the elliptic curves satisfying the condition of Theorem 11.31 is an infinite 
family. Finally in Section [7] we compute the bounds of z(P), which are used in Section 

m 



2. Preliminaries 

The standard symbols Q, K, C and Z will denote respectively the set of rational, 
real and complex numbers and the rational integers. We denote the discrete valuation 
on Z at the prime p by v p (-). We denote the set of all places of a number field K by 
M K . 

Throughout this paper, we assume that a, b G Z, a, b > 3, gcd(a, b) = 1 and 
m = a 6 + 166 6 . 
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As usual we write the Weierstrass equation for elliptic curves E over a number field 
K as 

(2.1) E : y 2 + a^xy + a 3 y = x 3 + a 2 x 2 + a 4 x + a 6 (a 1; a 2 , a 3 , a 4 , a 6 G X). 

Since the characteristic of K is not equal to 2, by completing the square of the left- 
hand side we have 

(2.2) (2y + a x x + a 3 ) 2 = 4x 3 + b 2 x 2 + 2o 4 :r + 6 6 , 
where 

b 2 = a\ + 4a 2 , o 4 = 2a 4 + 0103, &6 = a 3 + 4ae, 
08 = a t ae + 4:a2(iQ — aia 3 a 4 + a 2 a 3 — a 4 . 
Further, by putting 

c 4 = b\ - 246 4 , c 6 = -6^6 8 + 366 2 6 4 - 2 1 66 6 , 

we have 

{108(2?/ + a x x + a 3 )} 2 = (36x + 36 2 ) 3 - 27c 4 (36x + 36 2 ) - 54c 6 . 
We also define the discriminant of E as 
(2.4) A = -& 2 & 8 - 8b 3 4 - 276 2 + 96 2 6 4 6 6 . 

Using the form (12.31) . we can write 

(2 V r(2P) - X " - - 2h&X - b8 

[ ] 1 j " 4x 3 + 6 2 x 2 + 26 4 x + 6 6 

for P = (x, y) G E. 

Next we define the canonical height, which is a powerful tool to consider the 
arithmetic of elliptic curves. Let E be an elliptic curve over Q and P = (x, y) G 
E(Q). If x = n/d and gcd(n, d) = 1, we define the naive height of P by h(P) = 
max{log \n\, log \d\} ([12], p. 202]) and the canonical height of P by 

Hp) - iim ^ 

m P- 248]). 

Remark 2.6. In our definition the value of h is twice of those in [12], [1] and [TT] . 

The canonical height has the following properties. 

• h(P) = if and only if P is a torsion point. 

• h(kP) = k 2 h(P) for all P G £(Q) and all fceZ. 

• /i is a quadratic form on E. 

For details see also (T2J Chapter VIII Section 9]. 

Our computations of the canonical height is done by using the local height. We 
recall the existence of the local height function as follows. 

Theorem 2.7. (Neron, Tate, [HJ p. 341]) Let K be a number field, v a place and 
K v its completion respect to an absolute value \ ■ \ v . Let E be the elliptic curve over 
K given by ( \2.1\) . Then there exists a unique function X v : E(K V ) \ — > M. which has 
the following three properties. 
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(1) For all P E E(K V ) with 2P ^ O, 

X V (2P) = 4X V (P) - 2 log \2y(P) + a lX {P) + a 3 \ v . 

(2) The limit limp^,.o{X v (P) — log \x(P)\ v ) exists. 

v-adic 

(3) \ v is bounded on any v-adic open subset of E(K V ) disjoint from O. 

The function X v above is called the local height function. If we have to specify 
the elliptic curve, we may use the notation such as Xe, v - The canonical height can 
be decomposed as the sum of local heights. The sum of the local heights for all 
archimedean (resp. non-archimedean) places is called the archimedean (resp. non- 
archimedean) part of the canonical height and denoted by hf(P) (resp. /ioo(P))- We 
only consider the case K = Q and in this situation, 

(2.8) h(P)=h f (P)+h 00 (P)= A P (P) + A oc (P). 

p:prime 

Let d E K and 

7-1/ / i\2 i III, II I /\3 i // i\2 , ill l 

h : [y ) + a l xy +a 3 y = {x ) + a 2 [x ) + a 4 x + a 6 
the elliptic curve obtained by making the substitution 

(2.9) x' = x + d, y' = y 
in ([HI]). Then 

CL\ — &\) &2 ~ &2 — 3c?, CI3 = CI3 — dCL\^ 

(2.10) 2 „ „ 
a 4 = a 4 — 2da 2 + 3d , a$ = a§ — da^ + d a 2 — d . 

Now let P E E(K V ) and P' = (x(P) + d,y(P)) E E'(K V ). It is clear that the map 
E(K V ) 9?4 P 1 E E'(K V ) is a group isomorphism. 

Lemma 2.11. In the situation above, we have \e,v{P) — Ae'^P')- 

Proof. To see this, it is sufficient to show that the function / : E'(K V ) — > M. defined 
by f(P') = Xe,v(P) satisfies the three properties of A\, in Theorem 12.71 
The property (1) follows from the equality 

2y' + ai'x' + a 3 = 2y + a\{x + d) + a 3 — da\ = 2y + a\x + a 3 . 

For the property (2), we have 

lim {f(P') - log \x'(P%} = lim {Xe,v(P) ~ log \x(P) + d\ v } 

P'^tO' P^O 
?;-adic ti-adic 



lim ^A^(P)-log|x(P)|„-log 



u-adic 



lim lX E>v (P)-log\x(P)\ v -\og 

u-adic 

lim{A^(P)-log|x(P)|,}. 

P— 

u-adic 



x{P) + d 



x(P) 
d 



1 + 



x(P) 



The property (3) is clearly satisfied. □ 
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3. Computing the canonical height 

Let E a £ be the elliptic curve (II .ip and Pi, P 2 , P3 the rational points on E a ^ defined 
in ' 

Proposition 3.1. If ab is odd, v 3 (b) = 1 and m is square-free, then the canonical 
heights of the points P\, Pi, P3 have the following bounds 

1 - 1 

- log m - 0.7441 < h{Pi) < - log m + 0.5409, 

-logm - 0.7579 < h(P 2 ) < - logm + 1.0515, 
3 3 

- log m - 0.5113 < h(P 3 ) < - log m + 0.5665. 
3 3 

Proof of Proposition \3.1\ We use the decomposition (I2.8P to estimate the canonical 

height. We first estimate the archimedean part h 00 (P)=\ 0D (P) by using Tate's series 

with Silverman's shifting trick (|llj). 

Let E be the elliptic curve defined by (12. ip . For P £ E(M), we put 

t = t{P):=l/x(P), 

(3.2) z = z(P) := 1 - M 2 - 2b 6 t 3 - b 8 t 4 , 

w = W {P) := At + b 2 t 2 + 2b 4 t 3 + b 6 t 4 , 

where b 2 , 64, b^, bg are as in ( 12 .3[) . Note that we have x(2P) = z(P)/w(P). By the 
property of the local height (Theorem 12.71 (1)) we have 

A 00 (2P) = 4A oc (P) - 2 log \2y(P) + a x x{P) + a 3 \. 

Then using (12. 2p . we have 

A 00 (2P) - log |x(2P)| = 4A 00 (P) - 2 log |2j/(P) + fll x(P) + a 3 | - log |a;(2P)| 

= 4Aoo(P) - log |4x(P) 3 + 6 2 x(P) 2 + 2b 4 x{P) + 6 6 | 
-log| 2 (P)/ W (P)| 

= 4{A 00 (P)-log|a;(P)|}-log|^(P)|. 

Putting /i(P) := ^(P) - log \x(P)\, 

/i(2P) = 4 yU (P)-log|z(P)|. 
So if we ignore the convergence, we have 

n=0 

In fact, by Tate's theorem ( |11| Theorem 1.2]) if there is e > such that |x(P)| > e 
for all P £ E(R), then for any P £ E(M), log |z(2 n P)| is bounded independently of 
n and therefore 

1 00 

^(P) = log \x(P)\ + - 4" n log k(2"P)|. 

n=0 

For d £ Q and P e E(R), the point P' = (x(P) + d, y(P)) is on the curve 

(3.3) E' : (y') 2 + a^x'y' + a 3 'y' = (x') 3 + a 2 \x') 2 + a 4 V + a 6 ', 
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where 

ai = di, a<2 = d2 — 3d, a 3 = a 3 — da±, 
a/ = a 4 — 2da 2 + 3d 2 , a 6 ' = a 6 — da 4 + cPa 2 — d 3 
as we saw in ( 12 . 1 Oj) . We similarly put 

t' = t\pi) :=l/x\P'), 
(3.4) z' = z'(P') := 1 - b' 4 (t') 2 ~ 2b' 6 (t'f - K(t')\ 

w > = w '(p') : = At' + b' 2 {t') 2 + 2b' A {t'f + b' 6 (t') 4 , 
where b' 2 , b' A , b' 6 , b' 8 are the values obtained by replacing a±, . . . , a 6 by a[, . . . , a' 6 in 

The reason why we make this substitution is that we obtain the Weierstrass model 
to which we can apply Tate's theorem above. We call this the shifting trick following 
Silverman. 




X 



Figure 1. Figure 2. 

y 2 + a x xy + a 3 y = (y') 2 + a! x x'y' + a' 3 y' = 

x 3 + a 2 x 2 + a A x + a 6 (x') 3 + a' 2 (x') 2 + a' 4 x' + a' & 

Now let E = E af) . If P E E(R), then x(P) > -m 1 / 3 . So if we take d such that 
d > m 1 / 3 , then x'(P') = x(P) + d > — m 1 / 3 + d > 0. Therefore the assumption of 
Tate's result is satisfied and we have 

oo 

\e>,oo(P') = log\x'(P')\ + -^4-log|/(2"P')l- 

n=0 

On the other hand, by Lemma [2.111 X^r ^(P') = A^ j00 (P). So we have 

1 oo 

Aoo(P) = Xe,oo(P) = log\x'(P')\ + -^4- n log|^ , (2 I1 P / )l- 

n=0 

Let us compute Aqo^) by this formula, taking d = 2a 2 + 4b 2 . Then the condition 
d > m 1 / 3 is clearly satisfied. So we have 

(3.5) A 00 (P 2 )=log|a;'(^)| + i ^4""log|z'(2^)|. 

n=0 
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We put X = a/b and compute x'(P^), z'(P£, z'(2PQ, z'(4PQ. Since x(P 2 ) = 2ab, 
the following hold. 

• x'(P^) = 2ab + 2a 2 + Ab 2 (see flO} for the coordinate of P 2 ) 

• z '(p^) = (X 8 - 2X 7 + 2X 6 + 8X 5 + 2X 4 + 16X 3 + 16X 2 - 32X + 32)/(2X 8 + 
8X 7 + 28X 6 + 56X 5 + 98X 4 + 112X 3 + 112X 2 + 64X + 32) 

• z'(2P 2 ') = (X 32 + 4X 31 + 2X 30 - 32X 29 + 124X 28 - 64X 27 + 144X 26 + 864X 25 + 
288X 24 + 1344X 23 + 9408X 22 + 2688X 21 + 8256X 20 + 34176X 19 + 16512X 18 + 
59904X 17 + 237600X 16 + 119808X 15 + 132096X 14 + 546816X 13 + 264192X 12 + 
344064X 11 + 2408448X 10 + 6881 28X 9 + 589824X 8 + 3538944X 7 + 1 1 79648X 6 - 
2097152X 5 + 8126464X 4 - 4194304X 3 + 1048576X 2 + 4194304X + 2097152) 
/(2X 32 - 16X 31 + 64X 30 - 96X 29 - 96X 28 + 864X 27 - 1376X 26 + 256X 25 + 
6800X 24 - 13120X 23 + 16640X 22 + 35200X 21 - 70400X 20 + 141184X 19 + 
121472X 18 -217600X 17 +878624X 16 +123904X 15 -570368X 14 +2603008X 13 - 
1312768X 12 -753664X n +6422528X 10 -6127616X 9 +1884160X 8 +8126464X 7 - 
12845056X 6 +8912896X 5 +5767168X 4 -16777216X 3 +16777216X 2 -8388608X+ 
2097152) 

• z'{AP' 2 ) = (X 128 - 8X 127 +2X 126 +384X 125 +- • -+38685626227668133590597632) 
/(2X 128 + 32X 127 + 208X 126 + 448X 125 + • • • + 38685626227668133590597632) 

In the computation of z'iPQ, z'{2P'<^), z'^P!,), we used the software PARI/GP (Ver- 
sion 2.3.4) [9]. The commands 

allocatemem (64000000) ; 
E=ellinit( [0,0,0,0,m] ) ; 
Ed=ellchangecurve (E , [1 , -d , , 0] ) ; 

zd=l-Ed . b4* (x+d) ~-2-2*Ed . b6* (x+d) ~-3-Ed . b8* (x+d) "-4 ; 
zd_a=subst(zd,m,a~6+16*b~6) ; 
zd_b=subst(zd_a,d,2*a~2+4*b~2) ; 

zd_lP2p=subst (zd_b , x , 2*a*b) ; 

zd_ lP2=subst (subst (zd_ lP2p , a , X) , b , 1 ) 

zd_2P2p=subst(zd_b,x,ellpow(E, [2*a*b,a~3+4*b~3] ,2) [1] ) ; 
zd_2P2=subst (subst (zd_2P2p , a , X) , b , 1 ) 

zd_4P2p=subst(zd_b,x,ellpow(E, [2*a*b,a~3+4*b~3] ,4) [1] ) ; 
zd_4P2=subst (subst (zd_4P2p , a , X) , b , 1 ) 

compute them. 

Since x'{P!^f /m, z'(P£), z'{2P£), z'{AP£) are functions of X, by elementary calculus 
we can compute their maximum and minimum. So we can find the following bounds. 

^log(4m) <logx'(P 2 ') < ^log(57.2218701m), 

(3 g ) -0.6637015 <4~ x log z'(P^) < 0, 

-0.0433217 <4- 2 log^'(2P^) < 0.1396289, 
-0.0363430 <4- 3 log2'(4P^) < 0. 
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For example, to have the bounds of logx'(P 2 ') as above it suffices to show 

4 AW ( = W + v + V) < 57 . 2218701 . 

m V x + 16 / 

We type (the above codes are needed) 

f2p=(2*a*b+2*a~2+4*b~2)~3/(a~6+16*b~6) ; 

f2=subst(subst(f2p,a,X) ,b,l) ; 

df2=deriv(f2,X) ; 

df 2n=numerator (df 2) ; 

fac2=f actor (df 2n*l . ) 

to find the factorization of the numerator of ((2X 2 + 2X + 4) 3 /(X 6 + 16))'. We do not 
write the output here, but it turns out that the fourth root (X = 1.6484223 • • ■ ) is the 
only positive root. Then x'{P^f jm = 57.22187008 ■ • ■ by the following command. 

subst(f2,X,-subst(fac2[4, 1] ,X,0)) 

Since lim x _ >0 (2X 2 +2X+4)7(X 6 +16) = 4 and lim x ^ 0O (2X 2 +2X+4) 3 / / (X 6 +16) = 8, 
we have the bounds for logx'^P^) as above. 

To compute the bounds of z'^P^), z'(2P 2 '), z'{AP!^), we proceed similarly. Note that 
if only a, b are real numbers, d = 2a 2 + Ab 2 > m 1//3 is satisfied. Then log \z'(2 n P£)\ 
has a finite value by Tate's theorem. So the denominators of z'{P£), z'(2P^), z'(AP£) 
do not have real roots. For z'{P!^}, we type 

dzd_lP2=deriv(zd_lP2) ; 
dzd_lP2n=numerator(dzd_lP2) ; 
facdz_lP2=f actor (dzd_lP2n*l . ) 

to find the factorizations of the derivative of z^PQ. Then we substitute all the values 
of the positive roots to 4 _1 logz'^) by the commands 

4"-l*log(subst (zd_lP2 ,X, -subst (f acdz_lP2 [Di , 1] ,X , 0) ) ) 

where Di = 4 because only the fourth root is positive. Since limx-s.0 4 _1 log z'{P!^) = 
and lim x ^ 00 4- 1 log2'(^ > 2) = 4" 1 log(l/2) = -0.173286795 ••• , we have the bounds 
for A" 1 log z'{P^). 

For z'(2PQ and z'{AP^) we type 

dzd_2P2=deriv(zd_2P2) ; 
dzd_2P2n=numerator(dzd_2P2) ; 
facdz_2P2=f actor (dzd_2P2n*l . ) 

4"-2*log(subst (zd_2P2 , X , -subst (f acdz_2P2 [D 2 , 1] ,X , 0) ) ) 

dzd_4P2=deriv(zd_4P2) ; 
dzd_4P2n=numerator(dzd_4P2) ; 
facdz_4P2=f actor (dzd_4P2n*l . ) 

4"-3*log (subst (zd_4P2 , X , -subst (f acdz_4P2 [D 3 , 1] , X , 0) ) ) 
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where O2 = 6, 7, 8 and LI3 = 8, 9, 12, 13, 16 since they are all the indices which corre- 
spond to the positive roots in factorizations up to multiplicity. Since limx-s.o log z'i^Pi) 
= lim M logz'(4P^) = 0, lim x ^ 00 4- 2 logz'(2P^) = 4" 2 log(l/2) =-0.0433216 •■ ■ 
and Iim x _ >00 4" 3 log2; / (4P2) = 4" 3 log(l/2) =-0.0108304 ••• , by comparing the val- 
ues obtained we have the bounds for 4 -2 log z'(2P 2 ) and 4~ 3 \ogz'{4P! 2 ). 

For the estimate of the remaining terms z'{2 n P^) (n > 3) we use the following two 
lemmas, which we shall prove in Section [71 

Lemma 3.7. Let d = 2a 2 + 4b 2 or d = 3a 2 + 4b 2 . Then z'(P') < 120.531634 for any 

PeE„ b (R). 



Lemma 3.8. (1) If d = 2a 2 + 4b 2 , then 0.062326 < z'(P') for any P G E aM 

(2) Ifd = 3a 2 + 4b 2 , then 0.038068 < z'(P') for any P G E aib (R). 

Remark 3.9. In general there is Silverman's bound of z'(P') ([TTJ Lemma 4.1]), 
which gives a bound dependent on a, b. In our case we find that there is a bound of 
z'(P') independent ofa,b. 

We continue the proof of Proposition 13.11 Since (1/4) Y^=3^~ n = 1/192, we have 

1 1 °° 1 
(3.10) log(0.062326) < - ^ 4" n log z'(2 n P! 1 ) < log(120.531634). 

n=3 

By 1ET5J) . flap) and (jOljl . we have 

^logm- 0.295724 < A oc (P 2 ) < - log m + 1.513566. 

To compute the non-archimedean part hf{P<i), we use Lemma \3. 191 which is proved 
in the next subsection. Recall P2 = (2ab, a 3 + 4b 3 ). So a, (3,5 in Lemma [3.191 corre- 
spond to 2ab, a 3 + 4b 3 , 1 respectively. Therefore 

h f {P 2 ) = ~log2. 

Since h(P 2 ) = Aoo(-P2) + hf(P 2 ), we have 

-logm - 0.7579 < h(P 2 ) < -logm + 1.0515. 
3 3 

We can estimate h(P\), h(P$) similarly by taking d = 3a 2 + 46 2 , 2a 2 + 4b 2 respec- 
tively. □ 

Remark 3.11. The shifting width is not necessary to be 3a 2 + 46 2 ,2a 2 + 4b 2 . We 
choose the width which give good enough bounds. We do not have an idea to determine 
the width which give the best bound. 

3.1. Non-archimedean part. In this subsection we compute the non-archimedean 
part of the canonical height, which was required in the proof of Proposition 13.11 To 
do this, we use [TT| THEOREM 5.2]. But the Weierstrass equation of the elliptic 
curve to which we apply this theorem needs to be minimal at p to compute X p . Let 
n G Z be sixth power free and E the elliptic curve y 2 = x 3 + n. Then the Weierstrass 
equation of E is global minimal if and only if n ^ 16 (mod 64) ([2j Corollary 5.6.4]). 
Since in Theorem 11.31 we assume that a 6 + 166 6 is square-free, we may assume that 
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n ^ 16 (mod 64), and that y 2 = x 3 + n is global minimal. Further we may assume 
that V2{n) = 0, v^in) = 0. The reason is explained by the next lemma. 

Lemma 3.12. Let a,b G Z with gcd(a,b) = 1. 

(1) Ifn = a 6 + 16b 6 , then 

v 3 (n) = 0. 

(2) Ifn = a 6 + 16b 6 and n £ 16 (mod 64), then 

v 2 (n) = 0. 

Proof. (1) For any k G Z, we have (3k ± l) 2 = 9k 2 ± 6fc + 1 = 1 (mod 3). So since 
gcd(a, b) = 1, we have (a 6 , 6 6 ) = (1,0) , (1,1) or (0, 1) (mod 3). In all the cases n ^ 
(mod 3). 

(2) We assume that v 2 (n) ^ 0, and deduce a contradiction. Since V2(n) ^ 0, a = 
(mod 2). Since gcd(a, 6) = 1, b ^ (mod 2). Since 6 is odd, b 3 is odd and we can 
write b 3 = 2k + 1. Then we have n = a 6 + 16(2Jfe + l) 2 = 64A; 2 + 64fc + 16 = 16 
(mod 64), since a 6 = (mod 64). This is a contradiction. □ 

In the next two lemmas we compute X 2 (P), X 3 (P) assuming that v 2 (n) = 0, v 3 (n) = 
respectively. 

Lemma 3.13. Let n(^ 16 (mod 64)) be sixth power free and E the elliptic curve y 2 = 
x 3 + n over Q. Let P = (a/5 2 , (5/5 3 ) (a, 0,8 G Z, 5 > 0, gcd(a, 5) = gcd(0, 5) = 1) 
be a rational point on E. Assume v 2 (n) = 0. If v 2 (a) = 0, then X 2 (P) = 2v 2 (5) log 2. 
Ifv 2 (a) y£ 0, then A 2 (P) = -flog 2. 

Proof. Since n ^ 16 (mod 64), ?/ 2 = x 3 + n is global minimal. So we compute X 2 (P) 
following the algorithm (PH P-354, SUBROUTINE in THEOREM 5.2]). 

For the general Weierstrass equation (I2.ip and a point P on it, we put x := 
x(P), y := y(P). Further we define A, B, C, A for P as follows. 

A := v p (3x 2 + 2a 2 x + a 4 — aiy), B := v p (2y + a\X + 03), 

(3.14) C := v p (3x 4 + b 2 x 3 + 3b A x 2 + 3b 6 x + b 8 ), 
A:=X p (P)/\ogp. 

This is the same definition as in the algorithm but the value of A is twice of that in 
the algorithm. Recall that in our definition the value of the canonical height is twice 
of that in [IT]. 

For our elliptic curve, since a\ = a 2 = a% = 04 = 0, b 2 = b 4 = 6 8 = and b e = An, 
we have 

, f3a 2 \ (20\ f 3a(a 3 + An5 6 ) \ 

(3.15) A = v p ^— j , B = v p ^— J , C = v p ^ J . 

Note that C4 = (i.e. t>p(c 4 ) 7^ 0). This condition has an effect in the algorithm. 
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On this condition, by the algorithm we have 

2 max jo, ~v p (a/5 2 )^ if A < or B < 



(3.16) A 



IB 

— — if A, B > 0, C > 3B 

C 

-— if A B > 0, C < 3B 
Now we consider the case of p = 2. If 1)2(01) = 0, then A < and by (13.1 6p 

A 2 (P) = A log 2 = 2max jo,-^ 2 (a/5 2 )J - log 2 = 2v 2 (6) log 2. 

We assume that v 2 (a) 7^ 0. Then v 2 (5) = 0, since gcd(a, 5) = 1. So A, B > 0. 
Since P is on E, we have the equation n5 6 = j3 2 — a 3 . Since v 2 (n) = 0, t> 2 (/3) = 0. 
So B = v 2 (2j3) = 1 and C = v 2 (a) + v 2 (a 3 + 4n5 6 ) > 1 + 2 = 3. So C > 3B, and by 

A 2 (P) = A log 2 = -^log2 = -jjlog2. 

□ 

Lemma 3.17. We consider the situation of Lemma \3.13l Assume v 3 (n) = 0. If 
v 3 (/3) = 0, then A 3 (P) = 2v 3 (5)\og3. Ifv 3 (/3) ^ 0, then A 3 (P) = -|log3. 

Proof. We compute A 3 (P) following f)3.15p . (13 . 16[) for p = 3. 
If v 3 ((3) = 0, then B < and by fl3TT6|) 

A 3 (P) = A log 3 = 2 max jo, -h, 3 (a/5 2 )^ ■ log 3 = 2v 3 (<f) log 3. 

The last equality is as follows. If v 3 (S) = 0, then max {0, — ^v 3 (a/5 2 )} = 0. So 
max {0, — \v 3 (a/5 2 )} = v 3 (5). If v 3 (5) ^ 0, then since gcd(a,5) = l,i7 3 (ct) = 0. So 
max {0, -\v 3 (a/5 2 )} = v 3 (5). 

We assume that v 3 ((3) ^ 0. Then v 3 (5) = 0, since gcd(/3,<5) = 1. So B = 
v 3 (2/3/5 3 ) = v 3 (f3) > and A = v 3 (3a 2 /5 A ) = v 3 (3a 2 ) > 0. Since P is on E, 
n5 6 = f3 2 — a 3 . Since v 3 (n) = 0, v 3 (a) = 0. Using the equality a 3 + 4n<5 6 = /3 2 + 3n<5 6 , 

C7 = W3 (3a) + v 3 (a 3 + 4n5 6 ) = v 3 (3a) + w 3 (/3 2 + 3n5 6 ) = 1 + 1 = 2. 

So we have 3P > C. By ( 13TTB]) 

A 3 (P) = A log 3 = log 3 = -^log3. 

□ 

Lemma 3.18. Lei n 6 Z ie square-free and E the elliptic curve y 2 = x 3 + n over Q. 
Let P = (a/5 2 ,/3/5 3 ) (a,/3,5 E Z, 5 > 0, gcd(a,<5) = gcd(/3,5) = 1) be a rational 
point on E. We assume that p ^ 2, 3 . T/ien A P (P) = 2t> p (<5) logp. 
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Proof. We compute X P (P) following (ETT5jl . f l3TTSj) . At first if v p (a) = or = 0, 

then since 5 is an integer, A < or B < 0. So 

A P (P) = Alogp = 2 max jo, -^ P (a/5 2 )| • \ogp = 2v p (S) log p. 

The last equality follows from the same reason as that in the proof of Lemma 13.171 

Next we assume that v p (a) > and v p (f3) > 0. Then v p (5) = because gcd(a, 5) = 
1. Since v p ((3 2 — a 3 ) > 1 and n5 6 = (3 2 — a 3 , we have v p (n5 6 ) > 1. But n is square-free, 
v p (n) = or 1. So this case does not happen. 

□ 

By the previous four lemmas, we have the following lemma. 

Lemma 3.19. Let n G Z be square-free and E the elliptic curve y 2 = x 3 + n over Q. 
Let P = (a/5 2 ,/3/5 3 ) (a,/3,5 G Z, 5 > 0, gcd(a,<5) = gcd(/3, 5) = 1) be a rational 
point on E. Then the non-archimedean part of the canonical height of P is as follows: 



where 



h f (P) = 2\og5 + \' 2 (P) + \' 3 (P), 

(v 2 (a) = 0) , 

A'(P) = { 2 

--log 2 Ma)^0), 



= 0) 



Proof. 



h f (P) = A 2 (P) + A 3 (P) + K(p) 



= A 2 (P) + A 3 (P) + logp 

= A 2 (P) - 2v 2 (5) log2 + A 3 (P) - 2v 3 (5) log 3 + 2 log JJ 

v 

= A 2 (P) - 2*; 2 (5) log 2 + A 3 (P) - 2*; 3 (5) log 3 + 2 log 5. 

By Lemmas 13.131 and 13.171 we see that A 2 (P) — 2v 2 {8) log 2 and A 3 (P) — 2v 3 (5) log 3 
are nothing but A 2 (P) and A 3 (P) respectively. 

□ 



4. Uniform lower bound 



In this section we compute a uniform lower bound of the canonical height (Propo- 
sition [O]), that is a lower bound of the canonical height independent of P G E(Q). 
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Proposition 4.1. Let n6Z and let E be the elliptic curve y 2 = x 3 + n over Q. Let 
P = (a/5 2 ,f3/5 3 ) (a, (3,5 e Z, 5 > 0, gcd(a,5) = gcd(/3,<5) = 1) be a rational point 
on E. We assume that n > 0. Then we have 



AooOP) > ^logn + - log 



+ 0.31494685. 



Proof. Recall that in our definition the value of the canonical height is twice of that 
in pp. By Algorithm 7.5.7 [TJ and (£2} 



(4.2) \ OB {P) = ^-lo i 

16 



1 / ^(P) 2 ^ 1 



n (n+l) 

where q = exp(2iriuj2/uJi), 9 = ^2^ =0 (—^) n q 5 sin {27r(2n + l)Re(^p)/o;i}, A is 
the discriminant of E, u\ and are periods of E such that uj\ > 0, Im(cj2) > and 
Re(u2/u)i) = —1/2 and zp is the elliptic logarithm of P. Recall that zp is the complex 
number in {tiUJ\ + t 2 uJ2 '■ < ti,t 2 < 1} such that p(zp) = x(P) and p'(zp) = 2y(P), 
where p is the Weierstrass p-function. 
Note that q is a real number since 

= exp ( — ixi — 27rlm ( — - ) ) = — exp ( — 27rlm 



LO\ J J \ \UJi 

By Definition 7.4.6 and Algorithm 7.4.7 in [1] 

2tt _i 2tt 

w i = -—, ; = = n 6 



AGM(2 v / 3n5, \/2V^-3ns) AGM(2 v / 3, v^v^- 3) ' 

where AGM(-, ■) is the arithmetic geometric mean. So if we let ou[, u' 2 be the pe- 
riods of the elliptic y 2 = x 3 + 1, then we have u\ = n~e x u[. It turns out that 
u[ = 4.206546315 • • • . This can be done by PARI/GP (Version 2.3.4) ([9]) as follows. 

El=ellinit( [0,0,0,0,1]) ; 
El . omega 

Similarly by [U Algorithm 7.4.7], we have 

1 i AGM(2 v / 3ni V 2 V3 + 



U2/0J1 



2 2 AGM(2 v / 3ni, ^2^3 - 3ns) 

1 i AGM(2 v / 3, v / 271T3) , . , 

= 1 ; = 

2 2 AGM(2 ^3 ? ^2^-3) 

and so it turns out that q = —0.163033534 ■• • by PARI/GP as follows(the above 
commands are needed). 

-exp (-2*Pi*imag (El . omega [2] /El . omega [1] ) ) 
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Substituting these values and A = — 432n 2 in (14. 2ft . we have 



> 



16 
1 

16 



log 
log 



432n 2 



q 



log 



2tt<5 6 



log |0| 



432n 5 



0.163033535 



+ 4 log 



4.206546315n"6/3 2 



2n5 e 



log 1 1.167385748 1 



1 , 1, 

— log n H — lo£ 

12 S 2 6 



0.3149468597 •• 



by the trivial bound \9\ < 1 + \q\ + \q\ z + |(/| 



<?l 10 + \q 



15 



121 



< 1 



9 



M 1 



=1.16738574713 



□ 



Proposition 4.3. Let n be a positive, square-free integer and E the elliptic curve 
y 2 = x 3 + n. If P is a rational, non-torsion point on E, then 

(4.4) h(P) > — logra - 0.147152. 

Proof. By Lemmas 13.131 13.171 13.191 and Proposition 14.11 we have 

h(P) = h f (P) + \ 00 (P) 



> 2 log 5 + A' 2 (P) + A' 3 (P) + — log n + - log 



0.31494685 



= ~ log 5 + A' 2 (P) + { A 3 (P) + i log |0| J + ^ logn + 0.31494685 

> — log n - - log 2 + 0.31494685 = — log n - 0. 1471512 • • • , 
12 3 12 

since 5 G Z and A' 3 (P) + § log > 0. 

5. Estimate of the lattice index 



□ 



Let E be an elliptic curve of rank r(> 2) defined over a number filed K. Let 
Qi, Q2, Q s (s < r) be independent points in E(K). Then there exist generators 
G 1 ,G 2 ,...,G r of the free part of E(K) such that Qx, Q 2 , Q s e ZGi+ZG 2 + - ■ -+ZG S 
by the elementary divisor theory. The index of the subgroup ZQi + ZQ2 H — • + ZQ S 
in ZGi + ZG 2 + • ■ ■ + Zl7 s is called the lattice index of {Qi, • ••> Qs}- We put 

(Qi, Qj) = \ (hQi + Qi) - HQi) - HQj] 

R(Qi, Qi, Qs) = det ((Qi, <2j))i<;j< s • 

It is known that the canonical height h is a positive definite quadratic form on 
E(K)/E(K) tors . When we identify E(K)/E(K) tms ~ ZGi + ZG 2 + • • • + ZC7 r as 
Z-modules, h is the quadratic form defined by the symmetric matrix ((Gi, Gj))i<i,j<r- 
Let /(x) = Y^Ji j=i fi,j x i x j be a positive definite symmetric quadratic form. Then 
it is known that there exists a constant 7 n called the Hermite constant such that 

inf f(m) < 7„ det(f, A 
mez<-\{0} v 7 _ w ,3J 
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For example, 

jl = 1, 7 2 = 4/3, 7 3 3 = 2, 74 4 = 4, . . . 

In this section we estimate the lattice index. For this we use the following theorem 
of Siksek. 

Theorem 5.1. ([10j Theorem 3.1]) Let E be an elliptic curve of rank r (> 2) defined 
over a number field K . Let Q 1 , Q 2 , Q s (s < r) be independent points in E(K) and 
v the lattice index of {Qi, Q 2 , ■ ■, Qs}- Suppose that A > is a constant such that any 
point P e E(K) of infinite order satisfies h(P) > A. Then 

v<R(QuQ2,:,Q s ) l/2 (%/\y /2 - 

Proposition 5.2. Assume that m = a 6 + 166 6 is square-free, ah is odd and the 

discrete valuation v 3 (b) equals 1. If m > 6.38 x 10 22 (this is true for either a > 6321 
or b > 3982), the lattice indices of {Pi,P 2 }, {P2,P^\, {P3,Pi} are less than 5. // 
m > 19088 (this is always true), the lattice indices of {P±, P2} , {P2,Pz\, {-P35-P1} are 
less than 7. 

Proof. In this situation Pi,P 2 ,Pz are independent by Proposition 16.71 in the next 
section. Let A = ^logm — 0.147152. Then h(P) > A for any non-torsion point 
P e E a<b (Q). Now by Theorem EH it suffices to show that R(P U P j ) 1 / 2 ( 72 /A) 2 / 2 
is less than 5 or 7, when m > 6.38 x 10 22 or m > 19088 respectively for i ^ j 
= 1, 2, 3). Since 

1 = h(P n \h(P„\ - 

we have 



R(P 2 , P 3 ) = h(P 2 )h(P 3 ) - - A [h{P 2 + P 3 ) - h(P 2 ) - h(P 3 

4 h(P 2 )h(P 3 ) - \ {h(P 2 + P 3 ) - h(P 2 ) - h(P 3 )\ 



(P(P 2 ,P 3 ) 1 / 2 ( 72 /A) 2 / 2 } 



< 



3 A 2 
A h{P 2 )h{P 3 ) 
3 A 



2 



4 (| log m + 1.0515) (| log m + 0.5665) 
< 3 (J log m- 0.147152) 2 ' 

The last inequality follows from Propositions 13.11 and 14.31 By elementary calculus we 
see that the last bound is less than 25 if m > 6.38 x 10 22 , less than 49 if m > 19088 
and decreasing if m > e 2 . 

Since the upper bound of h(Pi) given in Proposition 13. II is less than those of h(P 2 ) 
and h(P 3 ), the cases of {Pi, P 2 }, {P 3 , Pi} are clear. □ 

6. Independence of P 1 ,P 2 ,P 3 



In this section we show that in the situation of Proposition 15.21 Pi, P 2 , P 3 are 
independent and the lattice index of {Pj, Pj} (i 7^ j) is not divisible by 2, 3. 

Lemma 6.1. Let n G Z and let E be the elliptic curve y 2 = x 3 + n over Q and 
Q E E{Q) \ P(Q)tors- We write x(Q) = u/s 2 with gcd(w,s) = 1. Then Q £ 2E(Q) 
in either of the following cases: 
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(1) n is odd, u ^ (mod 8) and s is odd, 

(2) n = 1 (mod 9), u = 2 (mod 3) and s ^ (mod 3). 

Proof. We assume that there exists R = (w/t 2 , z/t 3 ) £ -E(Q) with gcd(w,t) = 1 
such that Q = 2R and deduce a contradiction. By (12. 5p or the following PARI/GP 
commands, 

En=ellinit( [0,0,0,0,n] ) ; 
ellpow(En, [w/t~2,z/t~3] ,2) [1] 

we have x(2R) = (9w A - 8wz 2 )/(Az 2 t 2 ) and so u/s 2 = (9w A - 8wz 2 )/(Az 2 t 2 ). On the 
other hand (z/t 3 ) 2 = (w/t 2 ) 3 + n since R is on E. Eliminating z, 

(6.2) s 2 w(w 3 - 8nt 6 ) = Aut 2 (w 3 + nt 6 ). 

(1) If n and s are odd, then u> is even by (16.21) . Further t is odd since gcd(u> , t) — 1. 
Then t>2(w(u> 3 — 8nt 6 )) > 5 (note that if V2{w) = 1, w 3 — 8nt 6 = 8x even). So 
W2(4ut 2 (u; 3 + nt 6 )) > 5 and therefore V2(u) > 3. This is a contradiction since n ^ 
(mod 8). 

(2) Assume that n = 1 (mod 9), n = 2 (mod 3) and s ^ (mod 3). Note that if 
x ^ (mod 3), then = 1 (mod 9) (so modulo 3 also). 

Assume w = (mod 3). Then t ^ (mod 3) since gcd(w,t) = 1. So the left hand 
side of (16. 2p = (mod 3) and the right hand side of (16.21) ^ (mod 3). This is a 
contradiction. 

Assume w = 1 (mod 3). If £ = (mod 3), then the left hand side of (16.21) = 1 
(mod 3) and the right hand side of (16. 2p = (mod 3). This is a contradiction. 

If t ^ (mod 3), then the left hand side of (16.21) = 2 (mod 3) and the right hand 
side of (16.21) = 1 (mod 3). This is a contradiction. 

Assume w = — 1 (mod 3). If t = (mod 3), then the left hand side of (16.21) ^ 
(mod 3) and the right hand side of (16. 2p = (mod 3). This is a contradiction. 

Note that w 3 = —1 (mod 9). 

If t ^ (mod 3), then w 3 — 8nt 6 = (mod 9) and w 3 +nt 6 = (mod 9). So we can 
write w 3 -8nt 6 = 9Wi, w 3 + nt 6 = 9W 2 . Then by (Q we have s 2 w;-91^i = 4ut 2 -9W 2 
(mod 27). So s 2 wlVi = 4ut 2 iy 2 (mod 3). Therefore -W l = -W 2 (mod 3). On the 
other hand 9W 2 - 9Wi = 9nt 6 and so W 2 - Wi = nt s ^ (mod 3). This is a 
contradiction. □ 

Remark 6.3. Assume that we can write x(Q) = u/s 2 = u'/s' 2 (u',s' £ Z and not 
necessarily gcd(w', s 1 ) = 1). So u\u' and s\s' . Then if m' ^ (mod 8), u ^ (mod 8). 
If s' is odd, s is odd. If s' ^ (mod 3), s ^ (mod 3). If w' = 2 (mod 3), m = 2 
(mod 3) since u' = (s'/s) 2 u and s'/s ^ (mod 3). 

So it is not necessary to assume gcd(u, s) = 1 in Lemma [6. 1[ 

Lemma 6.4. Let n £ Z and /e£ 2£ 6e i/ie elliptic curve y 2 = x 3 + n over Q and 
Q £ £?(Q) \ J B(Q)tors- We write x(Q) = u/s 2 with gcd(w, s) = 1. T/ien Q ^ 3E(Q) 
m either of the following cases: 

(1) n is odd and u is even, 

(2) n = 1 (mod 9), u = 1 (mod 3) and f3(s) = 1. 

Proof. We assume that there exists i? = (w/t 2 , z/t 3 ) £ -E'(Q) with gcd(w,t) = 1 such 
that Q = 3R and deduce a contradiction. By the following PARI/GP commands 
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En=ellinit( [0,0,0,0,n] ) ; 
ellpow(En, [w/t"2,z/t~3] ,3) [1] 

we have x(3R) = (QAz 6 - U4w 3 z 4 + 81w 9 )/9t 2 w 2 (Az 2 - 3w 3 ) 2 and so u/s 2 = (QAz 6 - 
lAAw 3 z 4 + 81w 9 )/9t 2 w 2 {Az 2 - 3w 3 ) 2 . On the other hand (z/t 3 ) 2 = (w/t 2 ) 3 + n since 
R is on E. Eliminating z, 

(6.5) s 2 {(w 3 + Ant 6 ) 3 - 2 2 3Wi 6 } = 3 2 uw 2 t 2 (w 3 + Ant 6 ) 2 . 

(1) If u is even, then s is odd since gcd(w, s) = 1. Then since (w 3 + 4nt 6 ) 3 — 
2 2 3 3 nw 6 t 6 is even, it) must be even. So t is odd since gcd(w,t) = 1. Since n is odd, 
V2{w 3 + 4nt 6 ) = 2 and therefore u 2 (the left hand side of (16. 5p ) = 6. On the other 
hand v 2 (the right hand side of (16. 5j) ) > 7. 

(2) If v 3 (s) = 1, we can write s = 3s' (s' ^ (mod 3)). So we have 

(6.6) s' 2 {(w 3 + 4nt 6 ) 3 - 2 2 3Wt 6 } = uw 2 t 2 {w 3 + Ant 6 ) 2 . 

Now we show wt ^ (mod 3). Assume that wt = (mod 3). Then since the 
each side of (EBJ = (mod 3), we have (w 3 + Ant 6 ) 3 - 2 2 3 3 nw 6 t 6 = (mod 3). So 
w 3 + Ant 6 = (mod 3). But this does not happen since gcd(w,t) = 1 and n = 1 
(mod 9). So we see wt ^ (mod 3). 

Now if we assume that w = — 1 (mod 3), then w 3 + 4nt 6 = — 1 + At 6 = 3 (mod 9). 
So v 3 (w 3 + Ant 6 ) = 1. Then v 3 (the left hand side of flM])) > 3 
and t> 3 (the right hand side of (16. 6p ) = 2. This is a contradiction. 

If we assume that w = 1 (mod 3), then u> 3 + Ant 6 = —1 (mod 3). Then seeing 
(16. 6p modulo 3, we have u = — 1 (mod 3). This is a contradiction. 

□ 

Proposition 6.7. VFe assume that m = a 6 + 16o 6 is square-free, ah is odd and 
the discrete valuation v 3 (b) equals 1. Then P\, P2, P 3 , P\ + Pi, P2 + Pi, Pi + P 3; 
Pi + P2 + P 3 £ 2E a , b (Q) and P 1} P 2 , P 3 , P l ± P 2 , P 2 ± P 3 , P l ± P 3 , P x + P 2 ± P 3; 
Pi — P 2 ± P 3 ^ 3Ea t b(Q). In particular, Pi, P 2 , P 3 are independent and the lattice 
indices of {Pi, P 2 , P 3 }, {Pi,P 2 }, {P 2 ,P 3 }, {P 3 ,Pi} are not divisible by 2 nor 3. 

Proof. To ease the notation, we put P = P a ,6- We have 

x(P 1 ) = -a 2 , x(P 2 ) = 2a6, x(P 3 ) = -2ab, 

2a(a 3 + a 2 b - 2ab 2 - Ab 3 ) 



x{Pi + P 2 ) 
x(P - P 2 ) 

x(p + p 3 ) 

- P 3 ) 



(a + 26) 2 
2(a 4 - 3a 3 b + 6a 2 b 2 - 8ab 2 + 86 4 ) 
a~ 2 

2 (a 4 + 3a 3 6 + 6a 2 b 2 + 8a6 3 + 8b 4 ) 
a~ 2 

2a(a 3 - a 2 b - 2ab 2 + Ab 3 ) 



26) 



2 



4o 4 a 4 
x(P 2 + P 3 ) = -5-, x(P 2 -P 3 ^ 



a 2 ' v z a/ (26) 2 ' 
2a(a 5 + 4a 4 6 + 8a 3 6 2 + 12a 2 6 3 + 14a6 4 + 8b 5 ) 
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x{P 1 - P 2 - P 3 ) 



2a(a 5 - Aa% + 8a 3 o 2 - 12a 2 b 3 + Uab 4 - 8b 5 ) 
(a 2 - 2ao + 26 2 ) 2 ' 



Note that m = a 6 + 16o 6 = 1 (mod 9) since v 3 (b) = 1 and gcd(a, b) = 1. As we saw in 
Remark I6.3[ we can use Lemma 16.11 without the assumption that the x-coordinate is 
an irreducible fraction. Note that m in this proposition corresponds to n in Lemma 

o 

We see that Pi + P 2 £ 2E(Q) by Lemma 0(2) since 2a(a 3 + a 2 b - 2ab 2 - 4b 3 ) = 
2a A = 2 (mod 3) and a + 2b = a ^ (mod 3). Similarly P x + P 3 £ 2E(Q) by Lemma 
IOf2). It is clear that Pi, P 2 , P 3 , P 2 + P 3 , P x + P 2 + P 3 <£ 2E(Q) by Lemma 0(1). 

If there is a rational point R such that P\ = 3P, then h(Pi) = 9h(R). But by 
Proposition WM we have 9h(R) > 9(^logm - 0.147152) > |logm + 0.5409 > h(Pi) 
for m > 88, which is a contradiction. So Pi ^ 3P(Q). 

Since a 4 /(2o) 2 is an irreducible fraction, by Lemma [6.4( 2) P 2 — P 3 ^ 3P(Q). By 
computations we have 



x(2P 1 - 2P 2 - 
= a(-61446 17 



h 34816ao 16 - 101376a 2 o 15 + 204544a 3 6 14 
318248a 8 o 9 + 217216a 9 6 8 



320128a 4 o 13 - 
• 128160a 10 o 7 



439840a 6 o n + 403168a 7 o 10 

28152a 12 o 5 + 10200a 13 o 4 - 3006a 14 o 3 + 684a 15 o 2 - 108a 16 6 + 9a 17 ) /b 2 (2b 
40ab 5 + 56a 2 o 4 - 46a 3 6 3 + 28a 4 o 2 - 12a 5 b + 3a 6 ) 2 . 



409472a 5 o 12 
- 65072a n 6 6 
2 (166 6 



We denote the numerator by U' and the denominator by S' 2 . Further we write 
U' / S' 2 = U I S 2 as an irreducible fraction since it is an x-coordinate of an elliptic 
curve. Since v 3 (9a 17 ) = 2 and the orders of other terms of U' is greater than 2, 
v 3 {U') = 2. In 5", v 3 (b 2 ) = 2,i; 3 (3a 6 ) = 1 and other factors are not divisible by 
3. So v 3 (S' 2 ) = 4. Therefore, v 3 (S) = 1 and U" := U'/9, S" := S f /9 are integers. 
Clearly U" / S" 2 = U/S 2 . Since U" = a 18 = 1 (mod 3), U = 1 by the same argu- 
ment as in Remark 16 .31 So 2Pi — 2P 2 — P 3 G" 3P(Q) by Lemma [6.4( 2). Therefore 
Pi - P 2 + P 3 =-(2Pi - 2P 2 - P 3 ) + 3(P - P 2 ) £ 3P(Q). We have 



x(2P 1 + 2P 2 
- (40966 18 J 



hP 3 ) 

24576a6 17 



71680a 2 6 16 + 135680a 3 6 15 

.11 1 



188160a 4 6 14 



204800a 5 6 13 + 

181632a 6 6 12 + 133536a 7 6 n + 83488a 8 6 10 + 48472a 9 6 9 + 30720a 10 6 8 + 22464a n 6 7 + 
16496a 12 6 6 + 105 84a 13 6 5 +5496a 14 6 4 + 2 1 78a 15 6 3 +612a 16 6 2 + 108a 17 6+9a 18 )/a 2 6 2 (486 6 + 
128a6 5 + 156a 2 6 4 + 114a 3 o 3 + 56a 4 6 2 + 18a 5 6 + 3a 6 ) 2 



and by the same argument as above, we have 2Pi + 2P 2 + P 3 G" 3P(Q) by Lemma 
IBT4T2). Therefore Pi + P 2 - P 3 =-(2P 1 + 2P 2 + P 3 ) + 3(Pi + P 2 ) £ 3P(Q). We see 
that P 2 , P 3 , Pi ± P 2 , P 2 + P 3 , Pi ± P 3 , Pi + P 2 + P 3 , P - P 2 - P 3 G* 3E(Q) by Lemma 
16.4( 1). since the denominators of the x-coordinates of them are all odd. 

Next we prove the latter assertion of the proposition. By the elementary divisor 
theory there are generators Gi,...,G r G E(Q) and M G M 3 (Z) such that 





Pi " 




' G 1 ' 


(6.8) 


P 2 


= M 


G 2 




P 3 




G 3 
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Note that the lattice index of {Pi, P 2 , P 3 } equals | det M\. Let p be a rational prime. 
We have 



(mod pE(Q)), 



where M is the image of M in M 3 (Z/pZ). We assume that there exists A £ 
GL 3 (Z/pZ) such that AM has the row [0 0] and deduce a contradiction. Since 
we may assume that the first row is [0 0], by the left multiplication of A on (16. 9p 
we have 





Pi 




' G, ' 


(6.9) 


P2 


= M 


G 2 




P3 




G 3 





" hP, + k 2 P 2 + k 3 P 3 ' 




" " 




" G 1 ' 


(6.10) 


* 




* * * 




G 2 




* 




* * * 




G 3 



(mod pE(Q)), 



where [k\ k 2 k 3 ] is the first row of A. But the former assertion of this proposition 
implies that k x P x + k 2 P 2 + k 3 P 3 £ pE{Q) (p = 2, 3) for any (A*, k 2 , k 3 ) £ (Z/pZ) 3 \ 
(0, 0, 0). This is a contradiction. Therefore det M is not congruent to modulo 2 or 
modulo 3. 

By the same argument as above, the cases of {Pi, P 2 }, {P 2 , P 3 }, {P 3 , Pi} follow. □ 

Remark 6.11. By the same reason as above, if we verify that Pi, P 2 , P 3 , Pi ± P 2 , 
P 2 ±P 3 , P 3 ±Pi, Pi±2P 2 , P 2 ±2P 3 , P 3 ±2P 1 5E a>b (Q), we can prove that the lattice 
indices of {P U P 2 }, {P 2 ,P 3 } ; {p3,Pi} are not divisible by 5. Note that P $ 5P(Q) 
amounts to kP £ 5P(Q) (k = ±1,±2). For 3 < a < 6321, 3 < b < 3982 we can 
verify this by the software Magma ([7]). We give the code for this in AppendixlAl 

Now we can finish the proof of our main theorem. 

Proof of Theorem[TE For a > 6321,6 > 3982 by Propositions E21 EH the lattice 
indices of {Pi,P 2 }, {P 2 ,P 3 }, {P 3 ,Pi} equal 1. For 5 < a < 6321,3 < b < 3982 
by Propositions 15. 2\ 16.71 and Remark 16.111 the lattice indices of {Pi, P 2 }, {P 2 , P 3 }, 
{P 3 , Pi} equal 1. This completes the proof of Theorem 11.31 

□ 

We prove that there are infinitely many (a, 6)'s which satisfy the condition of The- 
orem [T31 



Lemma 6.12. The set 

S :-- 



m 



a b + lQb b £ Z 



a, b £ Z, m : square -free 
v 2 {ab) = 0, v 3 {b) = 1 



is an infinite set. 
Proof. We put 



S := < m = 



= (2k + 3/) 6 + 16(6Jfc - 9/) 6 £ Z 



k,l £ Z, m : square-free 



For {2k + 3Z) 6 + 16(6/c — 9/) 6 being square-free it is necessary that v 3 {k) = v 2 {l) = 0. 
Hence So is a subset of S. From the Theorem of Greaves ([5j THEOREM]) we see 
that S is an infinite set, since (2x + 3y) 6 + 16(6x-9?/) 6 = 8503785y 6 -34009308x?/ 5 + 
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56691900x 2 ?/ 4 -50384160x 3 ?/ 3 +25196400x 4 y 2 -6717888x 5 y+746560x 6 is an irreducible 
polynomial over Z. This is verified by the command factor of the software Maple 
([E])- Therefore S is an infinite set. 

□ 

7. Uniform bounds of z'(P') 

We use the notation of (I3.3p . (13 .4p . In this section we prove Lemmas 13.71 and 
13.81 which were used in Proposition 13.11 to give bounds of z'(P') independent of 

P e E a>b (Q). 

Although the following computations can be done by manually except the numerical 
evaluations, we included Maple commands to carry out all the steps at end of this 
section. 

Proof of Lemma \3. 7\ Let x = x(P). In this case since E' is y 2 = (x — d) 3 + m, we 
have a[ = a' 3 = 0, a' 2 = —3d, a' 4 = 3d 2 , a' e = m — d 3 , b' 4 = 6d 2 , b' 6 = 4m — 4d 3 , b' 8 = 
3d 4 - 12dm. By flU) 

, , Qd 2 Am - Ad 3 3d 4 - 12dm x 4 + Adx 3 - 8mx + Adm 

Z ^ )= ~ (x + d) 2 ~ 2 (x + rf) 3 {x + d) 4 {x + d) 4 ' 

Since x 3 + m = y 2 > 0, x > — m 1 / 3 . Note that d > m 1 / 3 , since (2a 2 + 46 2 ) 3 — (a 6 + 
166 6 ) = 486 6 + 96a 2 b 4 + 48a 4 6 2 + 7a 6 > 0. 
If x > 

x 4 + Adx 3 — 8mx + Adm x 4 Adx 3 Adm 

< 



If x < 



(x + d) 4 ~ (x + d) 4 (x + d) 4 (x + d) 4 

< l + 4 + 4 = 9. 

x 4 + Adx 3 — 8mx + Adm x 3 (x + Ad) —8mx + Adm 
{x + d) 4 = {x + d) 4 + {x + d) 4 

4 

— 8mx + Adm 8m 3 + Adm 
(x + d) 4 < (-mUd) 4 ' 
Assume d = 2a 2 + 46 2 . Putting Y = (a/6) 2 yields 

+ Adm 8(F 3 + 16) ((F 3 + 16)UF + 2 



4 

8m3 



(_ m | + rf )4 ^(F3 + 16) |_ 2 y_4) 4 

We denote the right hand side by g2,±(Y). Then 

d 48 (Y 2 - 8) (^2 y (y 3 + 16)^+4 (y 3 + 16)^ + 3 Y 3 + 48) 

^ 2 , 4 (y) = — -5 • 

ar (y 3 + 16)M (ys + 16)3 _2y-4j 

Note that (Y 3 + 16) 1 — 2 y — 4 < 0, since d — m 1 / 3 > 0. So 5'2,4(^ / ') has a minimum 
at y = \/8 and a maximum at y = — \/8. Therefore 

92,^00 ^ max I lim ^2,4 (y), lim g2,4(Y) \ ■ 
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Since 02,4(0) = 120.53163357 ••• , lim^oo #2,4(Y) = 16, we have z'(P') = # 2 ,4(Y) < 
120.531634. 

The case d = 3a 2 + 4b 2 is similar and we have z'(P') < 120.531634. □ 

Proof of Lemma \3.& We use the notation at the beginning of the proof of Lemma 
13.71 Let x = x(P), u = x/d and uq = —m l ^ 3 /d. Then u > uq > —1, since x > 
_ m i/3 > _^ p u tting Y = (a/b) 2 with substitution d = 2a 2 + 4b 2 yields 

, , x 4 + 4dx 3 — 8mx + 4dm d 4 u 4 + 4d 4 u 3 — 8dmu + 4dm 
Z ^ )= (x + d) 4 = (du + d) 4 

2it 4 (Y 3 + QY 2 + Y2Y + 8) + 8w 3 (Y 3 + QY 2 + 12Y + 8) - 2u(Y 3 + 16) + Y 3 + 16 
~ 2(w + l) 4 (Y + 2) 3 " 

We denote the last function by f(u,Y). Computing the derivatives, we have 

df 3(2u- l)(Y 2 -8) 



( M + l)4(y + 2 )4 ' 

df (4u 2 Y 3 + uY 3 - Y 3 + 24u 2 Y 2 + 48u 2 Y + 32u 2 + 16u - 16) 
3 



BY 

df 
du 



12(Y 3 + 6Y 2 + 12Y 



(u+l) 5 (Y + 2) 3 

{u-ui){u-u 2 ) 
} {u + lf{Y + 2f 



where 



Ml 



"2 



\/17Y 6 + 96Y 5 + 192Y 4 + 416Y 3 + 1536Y 2 + 3072Y + 2304 + Y 3 + 16 

8Y 3 + 48 Y 2 + 96Y + 64 ' 
V17F 6 + 96Y 5 + 192Y 4 + 416Y 3 + 1536Y 2 + 3072Y + 2304 - Y 3 - 16 

8Y 3 + 48 Y 2 + 96Y + 64 ' 



Clearly u\ < < u 2 for any fixed Y > 0, and the graph f(u, Y) has the either of the 
two forms below. Therefore /(it, Y) > min{/(u , Y), f(u 2 , Y)}. 



Figure 3. 



Figure 4. 
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At first we consider f(uo, Y). Since d = 2a 2 + 4b 2 , 

d 4 u 4 + 4d 4 u 3 - 8dmu + 4dm 



f(u ,Y) = f(-ms/d,Y) 



9m 3 



(du + dy 
9(r 3 + 16)3 



(d-mly (2Y + 4 - (Y 3 + 16)s) 4 

The last inequality follows from elementary calculus. 
Next we consider f(u 2 ,Y). 

df(u 2 ,Y) _ df 
du 



=-m3 jd 

> 0.75725080 



dY 



du2 df 



dY 
dY 



df_ 
dY 



3(2u 2 - 1)(Y 2 - 8) 
' (m 2 + 1) 4 (F + 2) 4 ' 



Now 
2u 2 



1 



2V17F 6 + 96F 5 + 192F 4 + 416F 3 + 1536F 2 + 3072F + 2304 - 2Y 3 - 32 



8Y 3 + 48Y 2 + 96Y + 64 



2V17F 6 + 96F 5 + ■ • ■ - (10F 3 + 48F 2 + 96Y + 96) 
8Y 3 + 48F 2 + 9QY + 64 



<0, 



since 



(Wl7r 6 + 96F 5 + ■ ■ -V - (10F 3 + 48F 2 + 9QY + 96) 2 
= -32F 6 - 576F 5 - 3456F 4 - 9472F 3 - 12288F 2 - 6144F < 0. 

So/(u 2 ,Y) > / (u 2 (V8), y/8). Since / (u 2 (V8), V8) = 0.06232685 ■■• , f(u 2 ,Y) > 
0.062326. Therefore z'(P') = f(u,Y) > 0.062326. 

The case d = 3a 2 + 4b 2 is similar and we have z'(P') > 0.03806854 • • • . □ 

The commands for those computations by Maple are as follows. 
The commands for Lemma [3. 7[ 

(8*m~ (4/3) +4*d*m) / (x+d) "4 ; 
subs (x=-nT (1/3) ,°/„) ; 
subs(m=a~6+16*b~6,d=2*a~2+4*b~2,°/o) ; 
g24:=simplify(subs(a=X,b=l,°/ )) : 
g24:=f actor (simplify (subs (X=sqrt(Y) ,%))) ; 
dg24 : =f actor (simplify (diff (°/„ , Y) ) ) ; 
evalf (subs(Y=0,g24) ,16) ; 
value (Limit (g24, Y=inf inity) ) ; 



(8*m~ (4/3) +4*d*m) / (x+d) "4 ; 
subs (x=-m~ (1/3) ,'/.) ; 
subs(m=a~6+16*b~6,d=3*a~2+4*b~2,°/ ) ; 
g34 : =s impl if y ( subs ( a=X , b=l , %) ) : 
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g34:=f actor (simplify (subs (X=sqrt(Y) //„))) ; 
dg34 : =f actor (simplify (diff (% , Y) )) ; 
evalf (subs(Y=0,g34) ,16) ; 
value (Limit (g34 , Y=inf inity) ) ; 

The commands for Lemma [3. 81 

(1) d = 2a 2 + 4b 2 

(x~4+4*d*x~3-8*m*x+4*d*m) / ( x+d ) ~ 4 I 
subs (x=u*d,°/„) ; 

subs (m=a~6+16*b~6 ,d=2*a~2+4*b~2 , '/.) ; 
subs(a=X,b=l,%) : 
subs (X=sqrt (Y) //„) ; 
f24 : =f actor (simplify (%)) ; 
df24Y:=f actor (diff (f24,Y)) ; 
df24u:=f actor (diff (f24,u)) ; 
sol24 : =solve (numer (df 24u) ,u) ; 
f24a:=subs(u=sol24[l] ,f24) : 
evalf (subs (Y=sqrt (8) , f 24a) , 16) ; 

9*nT (4/3) / (d-m~ (1/3) ) "4 ; 

subs (m=a~6+16*b~6 ,d=2*a~2+4*b"2 , %) ; 

subs(a=X,b=l, /.) : 

ff 24 :=f actor (simplify (subs (X=sqrt(Y) ,%))) ; 

factor (diff (%,Y)) ; 

evalf (subs (Y=sqrt (8) , f f 24) , 16) ; 

(2) d = 3a 2 + 46 2 

(x~4+4*d*x"3-8*m*x+4*d*m) / (x+d) "4 ; 
subs (x=u*d,°/„) ; 

subs (m=a~6+16*b~6 , d=3*a~2+4*b~2 , %) ; 

subs(a=X,b=l, /.) : 

subs (X=sqrt (Y) //„) ; 

f 34 : =f actor (simplify (%) ) ; 

df 34Y : =f actor (diff (f 34 , Y) ) ; 

df 34u : =f actor (diff (f 34 , u) ) ; 

sol34:=solve (numer (df34u) ,u) ; 

f34a:=subs(u=sol34[l] ,f34) : 

evalf (subs (Y=sqrt( 12) ,f34a) ,16) ; 

9*nr (4/3) / (d-m~ (1/3) ) "4 ; 

subs (m=a~6+16*b~6 , d=3*a~2+4*b~2 , %) ; 

subs(a=X,b=l,°/ ) : 

ff34:=f actor (simplify (subs (X=sqrt(Y) ,%))) ; 

factor (diff (°/„,Y)); 

evalf (subs (Y=sqrt( 12) ,ff34) ,16) ; 
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Appendix A. 

This is the code for Magma used in Remark 16. Hi By this code we can verify that 
P 1: P 2 , P 3 , Pt ± P 2 , P 2 ± P 3 , P 3 ± P u P ± 2P 2 , P 2 ± 2P 3 , P 3 ± 2Px i 5P a , b (Q) for 
5 < a < 6321, 3 < b < 3982 in the situation of Theorem 11.31 

for a in [5. .6321] do 

for b in [3. .3982] do 

if IsOdd(a*b) then 

if Valuation (b, 3) eq 1 then 

if Gcd(a,b) eq 1 then 

m:=a~6+16*b~6; 

if m It 6.381*10~22 then 

if IsSquaref ree(m) then 

E:=EllipticCurve( [0,0,0,0,m] ) ; 

P1:=E! [-a~2,4*b-3,l] ; 

P2:=E! [2*a*b,a~3+4*b~3,l] ; 

P3:=E! [-2*a*b,a"3-4*b"3,l] ; 

P4:=P1+P2; 

P5:=P2+P3; 

P6:=P3+P1; 

P7:=P1-P2; 

P8:=P2-P3; 

P9:=P3-P1; 

P0a:=Pl+2*P2; 

P0b:=Pl-2*P2; 

P0c:=P2+2*P3; 

P0d:=P2-2*P3; 

P0e:=P3+2*Pl; 

POf :=P3-2*P1; 

if 

{ 
< 

DivisionPoints (PI , 5) , DivisionPoints (P2 , 5) , DivisionPoints (P3 , 5) , 
DivisionPoints(P4,5) , DivisionPoints (P5, 5) , DivisionPoints (P6 , 5) , 
DivisionPoints (P7 , 5) , DivisionPoints (P8 , 5) , DivisionPoints (P9 , 5) , 
DivisionPoints (POa, 5) , DivisionPoints (POb, 5) , DivisionPoints(P0c,5) , 
DivisionPoints (POd, 5) , DivisionPoints (POe, 5) , DivisionPoints (POf ,5) 
> 

} ne 
{ 

<□, □, □, □, □, □, □, □, □, □> 

} then 

print a, b; 

end if ; 

end if ; 

end if ; 
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end if ; 
end if ; 
end if ; 
end for; 
end for; 
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